'''
Created on 11/05/2012

@author: nicolas
'''

import pylab as p

Fs = [20, 100, 101, 199, 201]       #frecuencia de muestreo
dTime = [1.0 / hz for hz in Fs]     #intervalo de muestreo

#defino el tiempo de muestreo
t_array = [p.arange(0, 1, time) for time in dTime]

#defino la funcion
def x(t):
    return p.cos(200 * p.pi * t)

#Muestreo
x_cos = [x(t) for t in t_array]

t_orig = p.arange(0, 1, 1.0 / 10000)
cos_orig = [x(t) for t in t_orig]


for i in range(len(x_cos)):
    p.figure(i)
    p.subplot(211)
    p.title('Original')
    p.plot(t_orig, cos_orig)
    p.plot(t_array[i], x_cos[i],'o')
    p.subplot(212)
    p.title(str(Fs[i]) + 'Hz')
    p.plot(t_array[i], x_cos[i])
    p.savefig('punto_1_cos_'+str(Fs[i])+'hz.png')

p.show()
